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Abstract 



We study the high energy behaviour of elastic scattering amplitudes within the leading 
logarithm approximation. In particular, we cast the amplitude in a form which allows us 
to study the internal dynamics of the BFKL Pomeron for general momentum transfer. 
We demonstrate that the momentum transfer acts as an effective infrared cut-off which 
ensures that the dominant contribution arises from short distance physics. 



1. Introduction 

We are interested in the study of elastic scattering in the Regge limit where s ^> — t. We 
might expect to be able to use the methods of perturbative QCD if the typical distances involved 
in the interaction are small. Whilst this may be true for some elastic scattering processes (for 
example, two highly virtual photons at large t) it is certainly not the case for others (for 
example, a pair of protons at t ~ 0). It is therefore necessary to check the self consistency of 
a perturbative calculation by investigating the typical distances involved in the interaction. If 
these turn out to be large then, even if the calculation is infrared finite, there is little justification 
for using a perturbative approach. In this paper, we consider the elastic scattering of generic 
colourless states via BFKL Pomeron exchange. We formulate the amplitude in such a way 
that we can readily identify the typical distances involved in the exchange. 

For elastic scattering at t — 0, the typical transverse momenta of the gluons which constitute 
the Pomeron are determined primarily by the sizes of the external particles. By studying the 
scattering of small size objects one therefore expects that perturbative QCD is valid. However, 
one must also take into account the diffusion properties characteristic of the BFKL exchange, 
namely as one moves away in rapidity from the external particles the width of the transverse 
momentum distribution broadens. This means that at sufficiently high centre of mass energies 
one is destined to pick up a large contribution from long distance effects. It is the main purpose 
of this paper to demonstrate that for t ^ 0, the situation changes rather dramatically. The 
scale \t\ effectively acts as an infrared cut-off, ensuring that the essential physics comes from 
the region where all the momenta are bigger than ~ \/— t. As a by-product we derive general 
formulae, which are useful in computing elastic scattering amplitudes in the Regge limit. 

2. The BFKL Pomeron 

We start by recalling the results of Lipatov p| for the solution of the BFKL equation for 
arbitrary momentum transfer, q 2 = —t. The Pomeron is described by the universal four-point 
function, F, depicted in Fig.l. In terms of this function, the amplitude for scattering the 
colourless states A and B elastically is 

A(s, t) = J d 2 k, d 2 k 2 $ A (k l5 q) F(y, k 1; k 2 , q) $ B (k 2 , q) . (1) 

The impact factors $a and determine the coupling of two gluons to the external states and 
we define them to contain no propagator factors (these are contained in F) . Rather than work 
in transverse momentum space, it is more convenient to work in the space of impact parameters. 
Accordingly we define 

Ffok^q) = / d 2 h, rfX d 2 h 2 d 2 b' 2 ( e -[k 1 .b 1+ (q-k 1 ).b' 1 -k 2 .b 2 -( q -k 2 ).m 

/fabbbi.ba,^)}. (2) 
For future notational convenience we denote this transformation 

F(y, ki, k 2 , q) = f {f(y, b 1; b' 1; b 2 , b' 2 )} . (3) 
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Figure 1: BFKL Pomeron exchange in elastic scattering. 



Lipatov determined that 

/ + °° V 2 — , N 

-oo ^ p + 1/4) 2 E " (bl ' ^ (b2 ' b ' 2) ^ 

where a$ = N c a s /ii and the function x(z/) is given by 

X (y) = 2^(1) - ^(1/2 + ^) - ^(1/2 - ^) • (5) 

Here ip(z) is the logarithmic derivative of the gamma function, T(z). The eigenf unctions, E u , 
are given by 

^i.bi) = [fe^r • (6) 



We have not included the contributions from non-zero conformal spin, since these vanish in the 
high energy limit. 

It is our aim to study the internal dynamics of the BFKL Pomeron. To do this we "break" 
the Pomeron into two pieces, as shown in Fig. 2. The procedure requires a little care, in order 
to account for the propagators which are present on the external legs. We utilize the following 
identity: 

F(y, k u k 2 , q)= J d 2 k F(y', k u k, q) F(y - y', k, k 2 , q) (7) 

where 

F(y, h, fc 2 , q) = hih - q)k* 2 (k 2 - q)* F(y, k x , k 2 , q). (8) 
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Figure 2: Elastic scattering amplitude as a convolution of the wavef unctions ipA and 

We have introduced complex numbers to represent the transverse momentum vectors (k = 
k x + iky). This definition means that F contains the part propagators l/(&i — <?)*, l/k 2 

and l/(k 2 — q) on its external legs. The wavef unctions, ipA,B, contain the BFKL dynamics and 
are given by 

F(y, k u k,q) 



MV, k, q) = / d 2 k, $ A (k 1; q) j^^f (9) 



and 



^(y, k, q) = / rf 2 k 2 $ B (k 2 , q) F t^ M) Q t (10) 



fc, fc 2 , q ) 

k* 2 (k 2 - q y 

The scattering amplitude is then given by 



A(s,t) = -^- Jd 2 k ^(y / ,k,q)^(y- y ',k,q). (11) 



(2tt)< 

In terms of Lipatov's solution in impact parameter space we have 



h(h - q) 

Similarly 

F(y, k u k 2 ,q) 



F u}t U -^ ) = f {^ d »0f(y> K b[, 62, b' 2 )} . (12) 



k :{k 2 _ q y = T {(^ K d bT )f(y, h, b[, b 2 , b> 2 )}. (13) 
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We can easily calculate d^dy f(y, &i, b[, b 2 , b' 2 ) using 



d b2 d b >E"*(b 2 ,W 2 ) = ^ ± ^ S-(b 2 ,b' 2 ). (14) 



We thus have 



fci, fc, g) _ . / f+oo ^ ^ g^X) E-(b 2 ,b 2 ) eHsx(l/) 1 (15) 



fci^-g) (z/ 2 + l/4) (& 2 -& 2 ) 5 

For our purposes it is beneficial to work in terms of the momentum transfer, q, and size, b. 
Transforming to this mixed representation, we define the new wavef unctions: 

<My, b,q) = /£ e- k - b V*(v, k, q) (16) 

and 

b, q) = / ^ e +4k - b ^(y, k, q) . (17) 
The scattering amplitude is now given by 

A(s, t) = | d 2 b * A (y', b, q) *%(y - y\ b, q). (18) 

By evaluating these wavefunctions we are able to compute the probability for finding the 
Pomeron with size b at rapidity y'. 

To this end, we shall proceed to compute the wavefunction, Using (|[|[15],[16|) we obtain 

Mv, b, q) = -^te J dv j^Jl) ^ (q) ^ (b ' q) (19) 

where 

V^(q) = j rf 2 k x rf 2 bx d 2 b; e -<(ki-bl+(q-ki).b' 1 ) $A ( ki; q ) b ,j ( 2 Q) 

and 

W„(b, q) = / ^ rf 2 b 2 rf 2 b' 2 e -*(k.b-k.b a -(q-k).bi) . (21) 

J 2tt (62 — 2 J 

All of the dependence on the impact parameter $^(k 1? q) is contained within the function 
V^ A (q). We proceed to calculate W„(b, q) and V^ A (q) separately. The k integral of ([2lD produces 
a delta function which fixes b = b 2 — b 2 . Introducing R = (b 2 + b' 2 )/2 we obtain 

WJ, b , 4 = %f*R e««™ ( (R _ b/2 ;; R + b/2)2 ) V2 "'" ■ (22) 
The integral over R can be performed using standard Feynman parametrization techniques: 

(27r) 2 |b| /q 2 V" f 1 j e " q ' b1 ' 



where K 2 i U is a modified Bessel function. The x integral can be performed exactly. In particular, 
we note that 



x 



/ dx K 2w (\ci\\b\Jx(l-x))=e-^ h / 2 / _=cos(-q.bVT^)tf aiv (-|q||b 

Jo ^ x (i _ x ) Jo VI - x 2 2 2 

(5 

and that the right hand side is tabulated in Gradshteyn & Ryzhik || (6.737(3)). We find 

2i7r 4 |b| / q 
r 2 (l/2-iz/) b 2 \ A 

where 

1 



^(b,q) = '"^i, 5 f^l V^A„(b,q) (2 ,) 



A,(b,q) 



sinh27rz/ 



t , v (q-b) 2 - q 2 |b| 2 - q.b ^ / ^/(q.b) 2 - q 2 |b| 2 + q.b 

6 "ii/ I . "ii/ 



(26) 



Equivalentlyf], we have 



2? 

A„(b, q) = — : — Im [J iu {bq* / A)J iv {b*q/ A)] . (27) 
smh27rv 



We next turn our attention to V^(q). This requires some choice for the impact factor 
$^(ki, q). We choose to write the impact factor in the form 

$ A (k, q) = J d 2 r f A (r) (e ik ' r/2 - e - lkr/2 ) (j^*' 2 - e -*(i- k W 2 ) . ( 2 8) 

This is a general choice, appropriate for scattering off colourless states. Assuming that /U(r) = 
/,i(|r|), the calculation of the function V^ A (q) can be computed using the method demonstrated 
in |J. The result is 



r(l/2 + ii/) \A J A/2-<oo 2th [7 V 16 

rri-z-zv) 



T(l/2 + «/2 - ii//2) T(l - z/2 - iv II) 

T(l-z + iu) 1 

T(l/2 + z/2 + iu/2) T(l - z/2 + iu/2) J ' 



(29) 



This is a particularly convenient way of expressing V A since it allows the limits 1/q 2 ^> r 2 and 
1/q 2 r 2 to be extracted with ease (by focusing on the poles lying closest to the contour). 
Alternatively, we could utilize (|24j) to write 



V u A (q) = 8*vr 5 J d 2 r f A (v) |r|e^ r K\ A„(r, q). (30) 



■""We note that this is consistent with the formula given in 
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This latter form does not assume /a(i") = /^(Irl). 



3. Inside the BFKL Pomeron 

We are now in a position to investigate the typical distances involved in the exchange. We 
can "look inside" the Pomeron and examine its size at some intermediate rapidity, y' . The 
probability density for having a Pomeron with size |b| at rapidity y' is proportional to the 
product of the wavefunctions, ^a{v', b, q)^s(y — y', b, q) (see (|18]) ) . We need to make some 
assumption regarding the nature of the external states in order to deduce the process dependent 
V,f-. We work with (|29|) under the assumption that the external state is characterised by a single 
scale, Q. In this case, we can write 



dr 



'q 2 r 2 ' 
~l6~ 



-f{r))=Q 



16Q 2 



h(z) 



(31) 



where h(z) is some dimensionless function of z which contains poles only to the left of the z 
plane contour, Re z = 1/2, otherwise the impact factor $^ is not defined. For photon and 
vector meson external states the following expressions for h(z) are needed: 



- f(r) = Q 2 K (Qr) 
-f(r) = Q 2 K 1 (Qr) 

By way of example we also note that 

ft \ Q -Q 2 r 2 

-f(r) = -e ^ 

-f{r) = rQ 3 e~ Q2r2 
-f(r) = 5(r 2 - 1/Q 2 ) 



Hz) 

h(z) 



2 2z - L T 2 (l/2 + z) 
2 2z ~ 1 T(l + z)T(z). 



(32) 



h(z) 

Hz) 
h(z) 



T[zl 
2 

r(i + ^) 

2 

1/2. 



(33) 



The existence of double poles induces additional logarithms in the final answer but does not 
affect any of our main conclusions. Subsequently, we will therefore assume that h(z) contains 
only single poles. We can now proceed to examine the wavefunction in the limits q 2 <C Q 2 and 
q 2 »Q 2 . 



The case when q 2 < 16 Q 2 

In this case the z-plane contour of (|29|) can be closed in the right half plane. We consequently 
pick up the poles in which come from the process independent part, i.e. at z — 1 ± iv. 
Evaluating the integral keeping only the nearest pole gives 



K(q) 



2(2vr) 5 Qr(l/2-iz/ 




h{\ + iv)T(-2iv) 
T(l/2 -iv)T(l + iv) 



(34) 
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The corrections to this expression are suppressed by powers of q 2 /g 2 . We have only the v 
integral to perform in order to obtain the wavefunction \I/ In order to proceed further we 
consider the two regions |q| > l/|b| and |q| < l/|b| separately. 

In the region where |q||b| < 1 we can use the small argument expansions of the Bessel 
functions, and for large y the dominant contribution to the v integral in flTSp will come from 
the region of small z/, hence we may write 

W,(b, q) » (£] W sin (v ln(16/(q 2 b 2 ))) . (35) 



2vb 2 \ 4 t 

Whilst for V v the corresponding small v limit of fl34|) is 



Qv V 4 
For the wavefunction, we thus obtain 



K A (q) - ( ^-Ysin (,ln(16Q 2 /q 2 )) . (36) 



*x(y,b,q) « 8^(1)^ / due asx(u)y sin (v ln(16/(q 2 6 2 ))) sin (z/ln(16g 2 /q 2 ) . (37) 



Expanding the BFKL eigenfunction about v = allows the integral to be performed. Using 
~ 4 In 2 - 14C(3)z/ 2 yields 



* A (y,b,q)=4^(l) 



Q6 2 {a 2 y)V 2 
Here a 2 = 14a s C(3). 



ln 2 (b 2 g 2 )\ /-ln 2 (256g 2 /(q 4 b 2 ))~ 



exp — exp 



Aa?y I \ Aa?y 



(38) 



In the region where |q||b| ^> 1 the function is again given by (|36]). For W u we take the 
large argument approximation to the Bessel functions. After some algebra we find that 



Inserting ( |3"B] ) and fl3Tf ) into our expression for (PD yields 

*Afe,b,q) « ^^(1 + e_iq ' b ) / ^ e 3 ^(">» z/sin (z/ln(16g 2 /q 2 )) (40) 
which evaluates (on expanding the eigenvalue) to 

V&A(y ' b ' q) = QHV (1 + e " tq - b )7^^ ln ( 16 ^ 2 /q 2 )exp(-ln 2 (16g 2 /q 2 )/4a 2 y). (41) 

Note that in the limit q^Owe only have the region described by ([38|) in which the second 
exponential term vanishes. This corresponds to the well known result for diffusion about the 
size of the external state, Q. In that case there is nothing to prevent diffusion into the large 



distance region where a perturbative calculation may be unreliable, (see, e.g. |]] for a more 
detailed study of diffusion at t = 0.) However, the effect of a finite q 2 is dramatic. Equation 
( |4"T| ) reveals that there is no diffusion inside the region |b||q| > 1. Moreover, the distribution in 
this region vanishes rapidly as b increases. The result is that diffusion to large sizes is effectively 
blocked beyond the scale l/|q| f7|. The momentum transfer acts as an effective infrared cut-off. 
In Fig. 3 we illustrate the diffusion for the case q = 0. This should be compared with Fig. 4 in 
which we show the equivalent distributions for the case q 2 = 1 GeV 2 . In generating these plots 
we chose a simple Gaussian form for the function f(r) which characterises the external impact 
factor, \f(r)\ = Q/re~ Q2r2 . As we noted earlier this corresponds to a function h(z) = T(z)/2. 

The case when q 2 > 16 Q 2 

For this case the z-plane contour integral must be closed in the left half plane. The only 
poles are those which arise from the process dependent factor h(z). In the limit q 2 ^> 16 Q 2 
the rightmost pole gives the leading power dependence. If this is a single pole at z = —z then 
we can take 

h( z ) » (42) 
z + z 

We can now evaluate Vj 4 : 



,Y{l/2-iu) fci 2 Y 1+lU ^ u , _ s ( 16Q 2 



K(q) = 2(2?r)5 1/n , - Mr Qh L (-z 



r(i/2 + ?v) v 4 J y V q 2 
r(i - z -w) 

T(l/2 + z Q /2 - iu/2) T(l - z /2 - w/2) 

T(l - z + iv) 
T(l/2 + z /2 + iu/2) T(l - z /2 + iu/2) ' 



(43) 



The corrections are suppressed by powers of Q 2 /q 2 . We can again separate the analysis into 
two regions consisting of large and small |q||b|. 

In the region where |q||b| < 1 the function W u (b,q) is again given by (p5|). For V v the 
corresponding small v limit of ([13]) is 




' " ] Qh L (z ). (44) 



q 2 



where we have further defined 



h ^-^(w^) 2 - (45) 

This form is only appropriate providing zq does not induce poles in the denominator, i.e. zq 
is not an odd integer. Although it would be straightforward to accommodate such values of 
zq (since they merely produce an additional factor of v 2 in the numerator) we shall ignore this 



S 



possibility in the following. This does not affect our main conclusions. For the wavefunction, 
we now obtain 

* A (y,b,q)*128h L (z )^^^) ° Jdvve^xMv sin (y ln(16/(q 2 b 2 ))) . (46) 

Expanding the BFKL eigenfunction about v = gives 

j / % \ aA r-1 , n Q|b| / 64Q 2 V° W1£!//2 , 2U /- ln 2 (16/(q 2 b 2 ))\ 
^(y,b,q) = 6V5FM«,)^i J p^ln(16/(q 2 b 2 ))exp ^ J . 

(47) 

In the region where |q||b| ^> I the function V„ is again given by (fy|) whilst W^(b, q) is given 
by fl39|). Inserting ( ^4|) and (|39|) into our expression for ^a (HD) yields 



* A (y, b, q) « 256^^(z ) j " (1 + e-*»' b ) / dv v 2 e a ^ (48) 

which evaluates (on expanding the eigenvalue) to 

O ~ /64O 2 V ■ ^ 

^(,,b,q) = 128 v ^^ L (, )^J (l + e-^)^,- (49) 

Notice that the wavefunction of fl4"T|) corresponds to diffusion about b 2 = 16/q 2 , i.e. the 
diffusion is no longer centred around the external scale, Q, determined by the impact factor. 
As before the momentum transfer q acts as an effective infrared cut-off. These properties can 
be seen in Fig. 5 

The wavefunction ty* B can be computed analogously. For the sake of illustration, we shall 
assume that the external states, A and B, are defined in terms of the scales Qa and Qb 
respectively and that in all other respects they are equal. In this case, is simply the complex 
conjugate of *&a with Qa replaced by Qb- Note that the complex factors in the denominators 
always combine to produce the factor 1/b 4 in the probability density (= ^a^*b)- We shall 
evaluate the probability distribution ^/a^*b f° r three separate cases. Firstly, when the external 
scales (Qa and Qb) are both larger than the momentum transfer; secondly, when one is large 
and one is small and finally when both are small. 

In Figures 6,7 and 8, we show the width of the 6-distribution as a function of y'. The width 
is determined by evaluating the value of lnb when b 2 ^A^B * s half of its maximum. Since 



A(s,t) 



J d(\nh 2 ) b 2 * A (y\ b, q)n(y - y', b, q) (50) 



it follows that for the dominance of short distance physics the "cigar" must remain inside 
the region of small |b|. Fig. 6 shows the case when both external scales are larger than the 
momentum transfer, q 2 . In this case the diffusion is similar to the q 2 = case, i.e. increasingly 
larger distances are probed as y' approaches y/2. However, movement to distances larger than 
~ l/|q| is blocked. This (almost) total exclusion from the region |b| > l/|q| remains even if 
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one or both of the external scales becomes small. In these cases the picture is quite different 
from the q 2 = case. Fig. 7 shows the diffusion when one external scale is large and the other 
is small. The case where both external scales are smaller than the momentum transfer, q 2 is 
shown in Fig. 8. 
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Fig. 3 
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Figure 3: The distribution of gluon size inside the BFKL Pomeron for \t\ = q 2 = and Q = 10 
GeV as described by ^ = b 2 /\h\. The plot shows the broadening of the distribution as the 
rapidity, y, is increased through the values y = 10, 15 and 20. Nothing prevents diffusion into 
the long distance region. The curves are from a numerical calculation of (|l^) using (|25|) and 
(|29|) . Each curve has been scaled by a factor e~^ oy '. 
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Fig. 4 
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Figure 4: The distribution of gluon size inside the BFKL Pomeron for \t\ = q 2 = 1 GeV 2 and 
Q = 10 GeV as described by ^ = b 2 /\h\. The plot shows the broadening of the distribution 
as the rapidity, y, is increased through the values y = 10, 15 and 20. The scale q acts as an 
effective infrared cut off preventing diffusion into the infrared region. The curves are from a 
numerical calculation of (|l^) using (p5|) and (p9|) . Each curve has been scaled by a factor e~ woy . 
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Figure 5: The distribution of gluon size inside the BFKL Pomeron for \t\ = q 2 = 1 GeV 2 and 
Q = 0.1 GeV as described by ^ = & 2 /|b|. The plot shows the broadening of the distribution 
as the rapidity, y, is increased through the values y = 10, 15 and 20. The scale q acts as an 
effective infrared cut off preventing diffusion into the infrared region. The curves are from a 
full numerical calculation of ( |19"D using ( p5|) and (29). Each curve has been scaled by a factor 
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Fig. 6 
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Figure 6: The width (as defined by half the maximum) of the distribution h 2 ^A^*B for the case 
Qa = Qb = 10 GeV and y = 20. The dotted lines represent the \t\ = case. The full lines 
represent the case \t\ = q 2 = 1 GeV 2 . 




Figure 7: The width (as defined by half the maximum) of the distribution b 2 ^^^^ for the case 
Q A = 10 GeV, Q B = 10~ 2 GeV and y = 20. The dotted lines represent the \t\ = case. The 
full lines represent the case \t\ = q 2 = 1 GeV 2 . 
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Figure 8: The width (as defined by half the maximum) of the distribution h 2 ^/^* B for the case 
Qa = Qb = 1CT 2 GeV and y = 20. The dotted lines represent the \t\ — case. The full lines 
represent the case \t\ = q 2 = 1 GeV 2 . 
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